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Abstract. 

A procedure avoiding any integration of the null geodesic equations is 
used to derive the direction of light propagation in a three-parameter family 
of static, spherically symmetric spacetimes within the post-post-Minkowskian 
approximation. Quasi-Cartesian isotropic coordinates adapted to the symmetries 
of spacetime are systematically used. It is found that the expression of the angle 
formed by two light rays as measured by a static observer staying at a given point 
is remarkably simple in these coordinates. The attention is mainly focused on the 
null geodesic paths that we call the quasi-Minkowskian light rays. The vector-like 
functions characterizing the direction of propagation of such light rays at their 
points of emission and reception are firstly obtained in the generic case where 
these points are both located at finite distances from the centre of symmetry. The 
direction of propagation of the quasi-Minkowskian light rays emitted at infinity 
is then straightforwardly deduced. An intrinsic definition of the gravitational 
deflection angle relative to a static observer located at a finite distance is proposed 
for these rays. The expression inferred from this definition extends the formula 
currently used in VLBI astrometry up to the second order in the gravitational 
constant G. 
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1. Introduction 

Testing gravity with projects like LATOR [I], ASTROD [2] and others will require a 
fully relativistic modelling of the direction of light propagation up tdj] terms of order 
G 2 , with G being the Newtonian gravitational constant (see, e.g., (3][3] and references 
therein). Such modelling appears to be relevant even for the Gaia mission [5], 
yet usually treated within the post-Minkowskian or post-Newtonian approximations 
( [EHS] and references therein). Indeed, an analytical solution for the direction of light 
propagation in the post-linear regime has been recently derived for a three-parameter 
family of static, spherically symmetric spacetimes with the view of explaining an 
apparent discrepancy between the standard post-Newtonian approach of Gaia and 
the numerical integration of the differential equations of the null geodesies [10] , 
Although being perfectly efficient, this solution presents the disadvantage of inferring 
the direction of propagation of light at its emission and reception points from the 
solution which corresponds to a ray emitted at infinity in a given direction. Our aim, 
in the present paper, is to furnish a new, more direct derivation based on a procedure 
developed in [11] and [12]. This derivation is natively adapted to the generic case 
where both the emitter and the receiver of the light ray are located at a finite distance 
from the source of the field. Furthermore, no integration of the null geodesic equations 
is needed. 

The paper is organized as follows. Section [2] lists the notations and conventions 
we use. Sections [3] and [3] yield the general formulae which characterize the direction 
of light propagation at the points of emission and reception in any static, spherically 
symmetric spacetime. The metric is assumed to be asymptotically flat and written 
for convenience in isotropic, quasi-Cartesian coordinates adapted to the symmetries of 
spacetime. In section [3] it is recalled how the direction of propagation and the impact 
parameter of a light ray can be derived from the time transfer function giving the 
travel time of a photon as a function of the spatial coordinates of the emitter and 
the receiver. In section 0] it is shown that the angle between two light rays passing 
through a given point as measured by an observer at rest at this point is equal to 
the Euclidean angle formed by the triples characterizing the directions of the rays 
when isotropic coordinates are used. In the subsequent sections it is assumed that the 
components of the metric can be expanded in power series of m/r, where m = GM/c 2 , 
with M being the mass of the central body and r the isotropic radial variable. The 
expansions are systematically stopped at order m 2 /r 2 . Furthermore, only the null 
geodesic paths for which the parametric equations can be expanded in power series of 
G are considered. These paths are called quasi- Minkowskian light rays in the present 
work. Section[S]is devoted to the calculation of the impact parameter and the direction 
of light propagation up to order G 2 in the generic case where both the emission and 
reception points are located at a finite distance from the centre of the gravitational 
source. The case of a ray emitted at infinity is treated in section [6] Section [7] yields 
an intrinsic expression of the gravitational bending of a light ray coming from infinity 
and observed by a static observer at a given point. This expression is the post-post- 
Minkowskian extension of the well-known formula currently used in VLBI astrometry. 
The classical expression of the light deflection in the post-linear regime is recovered 
in section [8] The conclusion is given in section [9] An appendix is devoted to the 
conservation law of the angular momentum in isotropic, quasi-Cartesian coordinates. 

| In the present work, "up to" always means "up to and including" . 



Direction of light propagation to order G 2 



3 



2. Notations and conventions 

We use the following notations and conventions. 

• Spacetime is assumed to be covered by a global quasi- Cartesian coordinate system 
(x 1 *) = (x°,x), where x° — ct, t being a time coordinate, and x — (x l ). 

• Greek indices run from to 3, and latin indices run from 1 to 3. 

• The signature adopted for the metric is H . 

• Any bold italic letter refers to an ordered triple. We distinguish the triples built 
with contravariant components of a 4-vector from the ones built with covariant 
components by systematically using the notations (a 1 , a 2 , a 3 ) — (a 1 ) = a and 
(6i, 62, 63) = (h) = b. All these triples are regarded as 3-vectors of the ordinary 
Euclidean space. 

• Given three triples a, b and c, c = a + b means that Cj = a' 1 + hi for any 
i = 1,2,3. Such a mixing of 'contravariant' and 'covariant' quantities cannot be 
ambiguous, since these equalities must be understood as equalities of triples, and 
not as equalities between true 4- vectors. 

• Given three triples a, b and c, a. b denotes the Euclidean scalar product a l b l and 
a.c denotes a l Ci, Einstein's convention on repeated indices being used in both 
cases. 

• \a\ denotes the formal Euclidean norm of the triple a: \a\ = (a. a) 1 / 2 . Similarly, 
|b| = (b.b) 1 / 2 . If \a\ = 1, a is conventionally called a unit (Euclidean) 3-vector. 

• a x b and axe arc the triples obtained by the usual rule giving the exterior 
product of two vectors of the Euclidean space. 

• The Euclidean angle between two triples a and b is denoted by S e (a, b). 
By convention, the determination of this angle is chosen so as to have < 
5 £ (a,b) < it. The angle S e (a, b) can be calculated by the relation S e (a, b) = 
arccos(a.6/|a|.|b|). The same notation and convention hold when a or b are 
replaced by 'covariant' triples c or d. Thus 5 £ (a,c) = arccos(a.c/|a|.|c|) and 
Se(c,d) — arccos(c.d/|c|.|d|). 

• Given two 4- vectors V and W,V-W stands for the Lorentzian scalar product of V 
and W: V ■ W — g(V, W) = g^V^W" , the quantities g^ v being the components 
of the metric tensor g. 

3. Direction of light propagation in a static, spherically symmetric 
spacetime 

Throughout this paper, spacetime is assumed to be a static, spherically symmetric 
Lorentzian manifold. For convenience, the metric is written in isotropic coordinates 
(x°,x) adapted to the symmetries of spacetime: 



where r = \x\ = \j8 i jx" l xi . The potentials A(r) and B(r) are supposed to satisfy the 
boundary conditions 




(1) 



lim A(r) = 1, 



lim B(r) = 1, 



(2) 
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which ensure that the metric is asymptotically flat when r — > oo. Moreover, we 
assume that there exists a value r^ > of the radial coordinate such that the metric 
is regular for any r > TV By analogy with the Schwarzschild spacetime, we suppose 
that r/i ~ m. The spatial region corresponding to r > r^ is denoted by T>h- 

Let x A and a; B be two points located in T>h- Generically, we consider a photon 
emitted at x A at a time t A = x° A /c and received at x B at a time t B = x B Jc, By 
convention, the instant of reception will be considered as fixed, whereas the instant 
of emission will depend on the trajectory of the photon. Owing to the static nature 
of the metric, any mention of t B will be systematically omitted and the null geodesic 
path followed by the photon will be denoted by T(x A , x B ) or simply T in the absence 
of ambiguity. A thorough characterization of the direction of propagation of the light 
at x A and x B may be carried out as follows. 

Let x a = a;" (C) be a system of parametric equations describing T(x A , x B ), £ being 
an arbitrarily chosen parameter. It is easily seen that the direction of propagation at 
any point x of T(x A ,x B ) is fully determined by the direction triple defined as 

l=(h/hM/hM/io) = {k/h), (3) 

where the quantities Iq and U are the covariant components of the vector tangent to 
the ray at x, that is the quantities given by l a = g a p(xr(())dXr(C)/d(- It is worth 
noting that the direction triple I is collinear to the Euclidean 3- vector tangent to the 
photon trajectory. Indeed a relation as follows 

T= 1 ^ (4) 

A(r)B(r) dx° 1 ; 

directly results from ([1} and ([3]). It is clear that I does not depend on the choice of 
the parameter along T(x A ,x B ). 

The expressions of I at points x A and x B will be denoted by I e {x A , x B \ T) and 
I r (x A , x B ; r), respectively. These triples may be determined, at least in principle, by 
integrating the null geodesic equations. However, the direction triples at x A and x B 
can also be straightforwardly derived from the relations [TT] 

T.(^.;r) = (ii) (5o , 



and 



£, (!r , a!jl :r, = (>) =-(c d%[x d A /«> ). m 



where Tr(x A ,x B ) is the so-called time transfer function associated to T, that is the 
expression yielding the travel time of the photon along T as a function of the spatial 
positions of the emitter and the receiver: 

t B - t A = T r {x A ,X B ). (6) 

Since the metric is static and spherically symmetric, 7~ r (x A , x B ) may be considered 
as a function of r A = \x A \, r B = \x B \ and the angle between x A and x B . So defining 
the unit vectors n A and n B by 

n A = — , n B = — , (7) 

r A r B 

and then setting 

H = n A .n B , (8) 
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T r {x Ai X B j = 7r \J*A j r B 7 A*) ■ 

As a consequence, equations (|5 a|) and (|56[) may be written in the form 



T e (:c A ,a; B ;r) = c*^n A + J"* 4 X Q AB x ra A , (9a) 

y / p\ d% \n A xn B \d% 

l r {x A ,x B ;T) = -c- — n B +c — Q AB x n B , (96) 

or B r B dn 

where Q AB is the Euclidean unit vector defined by 

= n A xn B 

\n A x n B \ 

It is well known that the angular momentum vector defined as 

L = -xxT (11) 

is conserved along any geodesic of a static, spherically symmetric spacetime (see 
appendix A). Let us denote by L r (x A ,x B ) the conserved angular momentum 
corresponding to the null geodesic path T(x A ,x B ). Replacing I by I e (x A , x B ; T) and 
taking (|9al) into account, it is easily S66I1 tllclt Tj^\X ^4 , X £ J may be written as 

Lr{ x at x b) — br{ x At x b) Q ab j (12) 

where b r {x A , x B ) is defined by 

d% 

b r (x A ,x B ) = -c—\n A X n B \. (13) 



Finally, using equation (fl3f , the system of equations (|9a| and (|96)) reads 

l e {x A ,x B ;T)=c—^-n A -Q AB X n A , (14a) 

ar A r A 

l r (a; A ,a;B;r) = -c— ^n B - — Q AB x n B . (146) 
c'^' B t B 

The main interest of introducing systematically the angular momentum vector lies 
in the fact that the constant of the motion 6 r is an intrinsic quantity, the magnitude 
of which may be regarded as the impact parameter of the ray when conditions ([2]) are 
met |13) . Indeed, for a null geodesic admitting an asymptote, the magnitude of L is 
such that \L\ = lim^-too \x x dx/cdt\ since I — > —(dx/cdt)^ when \x\ — > oo. So 

|6| = | — x x 1 1 is the Euclidean distance between the asymptote to the ray and the line 
parallel to this asymptote passing through the origin of coordinates O as measured by 
an inertial observer at rest at infinity. 

Equations (|14a[) and (|146[) appear to have a quite simple form. However, in order 
to discuss the problems raised by the bending of light, it will be more convenient to 
express I e and I r as linear combinations of the unit vectors N AB and P AB defined by 

N AB = = -n B - -n A (15) 

\x B x A \ \x B x A \ \x B — X A \ 

and 

P AB = N AB xQ AB , (16) 
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respectively. Then the triples I e and l r may be written as 



dT* b 
c—^(N AB .n A ) -\N AB x n A 

Or A r A 



N A 



97r &r 

— \N AB x n A | H (N AB .n A ) 

ur A r A 

c-q — (N AB .n B ) H |AT AB x n a | 

c^-\N AB x n B \ ~{N AB .n B ) 

or B r R 



(17a) 



N A 



(17b) 



Equations (fT"3|) . (|17a|) and (|17&P show that determining the direction triples l e 
and I r comes down to knowing the three quantities cdl^/dfj,, cdT T /dr A and cdl~ T /dr B 
associated to the null geodesic path covered by the photon. However, the problem 
is far from being solved, since two given points x A and x B can generally be linked 
by an infinite number of null geodesic paths. Furthermore, there exists in general 
no procedure enabling to determine the time transfer function corresponding to each 
null geodesic path joining x A and xj^. As far as we know, these problems can be 
successfully tackled only for the restricted class constituted by the null geodesic paths 
that we call below the quasi-Minkowskian light rays. This question will be developed 
at length in section [SJ after having supplemented the above general considerations by 
a remarkably simple expression for the angular distance between two light rays passing 
through a given point as measured by an observer at rest at this point. 



4. Angle between two light rays observed at the same point 

The only quantity having a practical physical meaning in astrometry is the angle 
between two light rays, say T and T', simultaneously arriving at a given point x B as 
measured by an observer passing through this point. The aim of the present section 
is to get a coordinate-independent expression of this angle in any spacetime endowed 
with a metric given by ([1]). For the sake of simplicity, we restrict our attention to the 
case where the observer is staying at rest at point x B . Such a static observer will be 
denoted by tSjJj]. The angle between the rays as measured by another observer passing 
through x B will be inferred using, e.g., the formulae given in |15j or |16) . 

It is well known that the direction of propagation of a light ray T as seen by 
an observer passing through x B is defined by the projection of a 4-vector l B tangent 
to r at point x B on the local rest space of this observer. Since the projection of a 
null vector is always a non-zero spacelike vector, the direction of propagation along T 
relative to S B is completely defined by the unit 4-vector (see, e.g., p"7TU9j ) 

n B = -^--U B , (18) 

U B ' "B 

§ This difficulty is facing us even in the Schwarzschild spacetime, in spite of the fact that the 
exact analytical solution to the differential equations of the null geodesies is known in terms of 
the Weierstrass elliptic function p for a long time (see, e.g., |14|). 

| The unit 4- velocity vector of S B is collinear to the timelike Killing vector at x B . So the property 
to be a static observer is an intrinsic one. 



Direction of light propagation to order G 2 



7 



where U B is the unit 4- velocity vector of S B , and U B ■ l B denotes the scalar product 
of U B and l B (see section [2]). Of course, the direction of T' arriving at x B following 
a tangent vector l' B will be defined by the 4-vector N' B obtained by substituting l' B 
for l B into (TT51) . For the sake of brevity, the angular separation between T and T' as 
measured by the static observer S B will be henceforth designated by S B (T, T') without 
explicitly mentioning the observer. The value of the angle 5 B (T,T') is fixed without 
ambiguity by the relation 

cos5 B (r,r') = -n B -n' B (19) 

supplemented by the condition 

< 5 B (T,T') < 7T- (20) 

The notation used for the angles as measured by a given observer must be carefully 
distinguished from the notation of the Euclidean angles between the direction triples 
defined in section [5J Surprisingly, however, there exists a very simple connection 
between the covariantly well-defined angles and the formal Euclidean angles. We have 
indeed a proposition as follows. 

Proposition 1 Let T and V be two light rays arriving at x B at the same instant. 
When the metric is written in isotropic coordinates adapted to the symmetries of 
spacetime, the angle 5 B (Y,T') between these rays as measured by a static observer 

staying at x B is equal to the Euclidean angle between the triples I r and l r defining at 
x B the directions of propagation along T and V , respectively: 

6 B (TX) = 5 £ {trXr) = arccos ( ^4r~] ■ (21) 

Proof of proposition [TJ Since the contravariant components of U B are U B = 
\j \J A(r B ) and U l B = 0, equation (|T5)) implies that the covariant components of N B 
a nd N' B are given by {N B ) = 0, (N B )i = ^/A{^) {T B )i and {N' B ) = Q,(N' B )i = 
\J A(r B ) (l' B )i, respectively. Consequently, equation (|19j) may be written as 

cos<5 B (r,r') = A{r B )B{r B )T_ r A r . (22) 

It may be straightforwardly deduced from the property of l a = g a pl@ to be a covariant 
null vector that 

|f P | - 1 (23) 

y/A(r B )B(r B ) 



for any null path T. Equation (|22|) is therefore equivalent to equation ([21]) when 
condition (|20| is taken into account. QED. 

The angular separation <5 B (r, T') is an intrinsic quantity defined as an angle 
effectively measurable by an observer who is himself defined in a way which does not 
depend on the choice of coordinates. So the equality of <5 B (r, T') with a purely formal 
Euclidean angle S £ (l r ,l r ) was not obvious. In fact this equality would be broken in 
any coordinate system differing from the isotropic one. The numerical coincidence 
yielded by proposition [1] therefore confirms the relevance of our choice of coordinates. 
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5. Post-post-Minkowskian approximation 

From now on we assume that in the exterior region T>h the functions A{r) and 6(r) -1 
admit analytic expansions in powers of m/r as follow 

. , , 2m m 2 
A(r) = 1 + 2/3— + • • • , (24a) 

1 „ m 3 m 2 . 

= l + 27 _ + -£— + •.., 24& 
£>(r) r 2 r z 

where /3 and 7 are the usual post-Newtonian parameters, and e is a post-post- 
Newtonian parameter ( J 9 = 7 = e= lin general relativity) . Since m is proportional to 
the gravitational constant G, equations (|24a[) and (|246[) may be considered as yielding 
a post-post-Minkowskian expansion of the metric. In agreement with this assumption, 
we may suppose that given an arbitrary reception point x B in T>h there exists a large 
region of T>h such that any point x A in this region is linked to x B by one and only 
one null geodesic path which is a perturbation in powers of G of a Minkowskian 
null geodesic. We may admit that such a configuration correspond to almost all the 
situations effectively encountered in astrometry. To be more precise we shall adopt 
the following definition. 

Definition 1 A null geodesic path joining two points x A and x B in the domain of 
analyticity T>h is said to be a quasi- Minkowskian light ray if this path lies in T>h and 
can be described by parametric equations which take the form 

00 

x° = x° B - X\x B - x A \ + G n X° (n) (x A ,x B ,\), (25) 

00 

x = x B - X(x B - x A ) + ^2G n X {n) (x A ,x B ,\), (26) 

n=l 

where X is the affine parameter varying on the range < A < 1, with the boundary 
conditions 

X? n) (x A ,x B ,0) = (27) 



X^(x A ,x B ,0) = X {n) {x A ,x B , 1) = (28) 
being met whatever n > 1 . 

It is shown in |12j that the time transfer function corresponding to a quasi- 
Minkowskian light ray is yielded as a series in powers of G by an iterative procedure 
involving only integrations along the segment defined by parametric equations x = 
x B —X(x B —x A ) 1 with < A < 1. This procedure implies the uniqueness of the solution. 
The path between x A and x B corresponding to definition [T] will be henceforth denoted 
by T s (a; A , x B ) or simply by rfl. For the sake of brevity, we shall drop any mention 
of path T s in the notations previously introduced for the time transfer function and 
the impact parameter. Thus T(x A ,x B ) and b shall henceforth stand for Tr s (x A ,x B ) 
and £> r3 (x A , x B ), respectively. Moreover, we shall systematically omit any mention of 
T s in the arguments of the direction triples. 



% The index s is adopted in order to recall that a quasi-Minkowskian null geodesic path may be 
considered as an almost straight line. 
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Since each term in the perturbation expansion of T(x Al x B ) is an integral 
involving the values of the metric tensor on the segment joining x A and x B , we shall 
assume that the condition 



\x B — X(x B — x A )\ > r/j when < A < 1 



(29) 



is fulfilled in order to guarantee the validity of the expressions we use. It is worthy of 
note that condition (|2"9"j) is differently ensured according to the signs of N AB .n A and 
N AB .n B : 

a) If N AB .n A > 0, ([29)1 is met if and only if r A > r/,; 

b) If N AB .n A < and N AB .n B < 0, (|2"9")l is met if and only if r B > r h ; 

c) If N AB .n A < and N AB .n B > 0, (JUJ) holds provided that 

r c > r h , (30) 

where r c is the so-called zeroth-order distance of closest approach, that is the Euclidean 
distance of the straight line passing through x A and x B from the origin O: 



r A r R \n A x n r 



(31) 



We shall not examine here if the results established below might be considered as 
valid under a less constraining condition than (|2"9"|) . 

It is shown in and |12) that the time transfer function corresponding to 
T s (x Al x B ) is given 

(1 + j)m 



T(x A , x B ) — 



X R X : 



■hi 



c 

X R X A 



+ r B + \x L 



J a + r B - \xt 
arccos(n A .n B ) 



(1+7) 2 
1 + n A .n B 



+ 0(G 3 ), (32) 



r A r B c 

where k is the dimensionless parameter defined by 

K=i[8(l+ 7 )-4/? + 3e]. (33) 
It may be noted that k = 15/4 in general relativity. 

The impact parameter of T s (x A ,x B ) is easy to determine. Substituting for 
T(x A ,x B ) from (|32|) into (TT3)) . noting that 

\x B - x A \ = (r\ + r\ - 2r A r B fi) 1 / 2 , 

and then using the relations 



r A n A x n F 



X R X 1 



= \N AB x n B \, 



\X fj X £ 



= \N AB xn A \ (34) 



and 



n A .n B - \N AB x n A \\N AB x n B \ = (N AB .n A )(N AB .n E 
a straightforward calculation leads to the following proposition. 



(35) 



+ The term of order G in (132 I t is the well-known Shapiro time delay (see, e.g., 1201 and references 
therein) . 
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Proposition 2 Let x A and x B be two points such that condition h29\) is fulfilled. The 
quantity b associated to a quasi- Minkowskian light ray T s (x A ,x B ) by equation U3\) is 
given by 



1 



(1 + 7)m \N AB x n A \ + \N AB x n E 



1 + n A .n B 

arccos(n A .n £ 



-(N AB .n A )(N AB .n E 



-(1+7) 



2 1 - {N AB .n A )(N AB .n E 



O 



l + n A .n B j 

where r c is the zeroth-order distance of closest approach defined by equation [37 



(36) 



Since the right-hand side of (1361) is a positive quantity when r c 3> m, we 
shall henceforth identify b with the impact parameter of the quasi-Minkowskian ray 

We are now in a position to determine the triples l e and l r corresponding to 
T s (x A ,x B ) by substituting for T(x A ,x B ) from (|32[) into (|17a|) and (|176[) . respectively. 
Using once again (|34|) . taking (|36|) into account and using the relation 

\n A x n B \ + \N AB x n B |A/" AB .n A - \N AB x n. 4 |AT AB .n B = 0, (37) 

a straightforward calculation leads to the following proposition. 

Proposition 3 Under the assumption of proposition the triples l e and l r 
corresponding to a quasi-Minkowskian light ray T s (x A ,x B ) admit the expansions in 
powers of mjr c given by 



l e {x A ,x B ) = —N A 



m\N AFt x n. 



m\N AB x rij 



1 + 7 + 



\N AR x n. 



>N A 



(1+7) 



-N AB .n A 



-(1+7) 2 



l + n A .n B r c 

N AB .n B - N AB .n A 
1 + n A .n R 



(1+7) 2 
1 + n A .n E 

arccos(n A .n s 
\n A x n B \ 

P AB +o(^\ (38a) 



-N AB .n r 



and 

l. r y^Ai 35b) = ~N AB 

m\N AB x n f 



m|AT 4R x n F 



1+7 + 



(1+7) 



m\N AB x n E 



n A x n F 



-N AB .n r 



+ (l + 7) 2 



1 + n A .n B r c 1 

N AB .n B - N AB .n A 
1 + n A .n R 



(1+7) 2 
1 + n A .n 

arccos(n A .n B ) 
\n A x n R \ 



>N A 



N AB .n A 



>Pab+0[ — ) , (386) 



respectively. 
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Owing to the intrinsic nature of the impact parameter, it is relevant to form the 
expansions of I e and l r in powers of m/b. We just have to replace l/r c into (|38ap and 
(|386p by its expansion inferred from (f3f))) . namely 



1 _ 1 

r c b 



(1 + jym \N AB x n A \ + \N AB x n E 
b l + n A .n B 



O 



rn 
1? 



(39) 



Using the relation 



[\N AB x n A \ + \N AB x n B \]±^ ^ + N AB .n A - N AB .n B = 0(40) 

l + n A .n B 

we obtain the fundamental expressions below. 

Proposition 4 The triples I e andl r yielded by proposition^ may be rewritten in the 
form 



Lei^Ai 2-b) — N AB 



m\N An x n, 



+ (1+7)2 |AUxn 



1 + n A .n 

m\N AR x n 



1+7 + 



>N A 



k\N ar x n. 



(1+7) 



|n A x n B 
1 + n A .n f 



Km 
~b~ 



l r (x A ,x B ) = —N A 



arccos(n A .ri B ) 
|n A x n B \ 

m\N AR x nJ I 



N AB .n B - N AB .n A 



O 



m 



(41 a) 



2 ^xn, 
1 + n A .n E 



m 

1 + 7+y 



' N A 



k\N ab x n F 



+ - 



m\N AB x n f 



(1+7) 



|n A x n B 
1 + n 4 .n F 



Km 
~b~ 



arccos(n A .n B ) 



N AB .n A - N AB .n £ 



O 



\n A x n 

where b is the impact parameter of the ray given by equation {5 



m 



(416) 



Equations (|41a[) and (|416[) have clearly a more simple form than equations p8a[) 
and (|386p , respectively. 

It may be noted that the Euclidean norm \l r \ is now given by 



\U = 1 



(l + 7 )m 1 



+ 2/c-(l+7) 2 -5-+° 



r B 2 

which implies that the factor l/|£ r |.|£ r | involved in (f^Tj) reads 

2(1 + 7)m 



1 

iLi.il 



= 1 - 



r B 



-2[«-2(l + 7 ) 2 ] 



m 



O 



m 



(42) 



(43) 
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6. Case of a ray emitted at infinity 

The notion of quasi-Minkowskian light ray can be extended as follows when the ray is 
emitted at infinity. 

Definition 2 Let N e be a unit 3-vector and x B a point in T>h- A null geodesic path 
is said to be a quasi-Minkowskian light ray emitted at infinity in the direction N e and 
received at point x B if this path entirely lies in 2\ and admits parametric equations 
having the form 

oo 

x Q =x° B +a-a B + J2 G n Xf n) (N e , x B ,a), (44) 

71=1 

OO 

x = x B + (a - a B )N e + G n X (n) {N e ,x B ,cr), (45) 

71 = 1 

where a is an affine parameter varying in the range — oo < a < o~ B , the boundary 
conditions 

X° n) (N e ,x B ,<r B ) = (46) 

and 

X {n) (N e ,x B ,a B ) = (47) 

being fulfilled. 

Such a ray will be denoted by r s (AT e , ai B ). Its impact parametei0 b is given by 
equation ([55)1 in which x A may be any arbitrarily chosen point on T s (N e , x B ). Indeed, 
it is easily seen that introducing the affine parameter A = (cr B — o~)/(o~ B — a A ) transforms 
the part of T s (N e ,x B ) joining x A and x B into a null geodesic path satisfying all the 
criteria set in definition [T] Consequently, b may be obtained as a function of 7V e and 
x B by taking the limit of the right-hand side of (f5o]l when r A — > oo, n A — > —N e and 
N AB — > N e . It will be useful to note that arccos(n_4 • n B ) — > ir — arccos(7V e • n B ) when 
x A is moved away along the ray towards infinity. It may be pointed out, however, 
that a condition like (|29|) has to be fulfilled by x A in order to guarantee the validity of 
(|36p . Unfortunately, it seems impossible to know if such a condition is effectively met 
by any point of T s (N e , x B ). So we shall henceforth content ourselves with assuming 
the validity of the asymptotic limit of condition (1291) , namely 



|(<7 — cr B )N e + x B \ > r/i when — oo < a < cr B . (48) 

Condition (|48l) means that the straight segment parallel to 7V e and ending at point 
x B is entirely inside the domain T>h- A proposition as follows can be formulated. 

Proposition 5 Let N e be a unit vector and x B a point in T>h chosen so that condition 
|^<S| ) is met. The impact parameter b of a quasi- Minkowkian light ray Y s {N e ,x B ) is 
given by 

, \, , (l + 7)m \N e xn B \ 

b = r c < 1 H — 

r c 1 - N e .n B 

* It is natural to assume that | X/tw=i ^™-^(n) C^ei x B , a)\n B — (N e .n B )N e ] | gcnerically admits a 
bounded limit when a — > — oo. This assumption means that the ray is asymptotic to a straight line 
parallel to iV e . 
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7r — arccos(iV e .n B ) 
1 + — V , ; A^ e -rt a 



(l + 7 ) ; 



. 1 + AT e .w £ 
1 - N P .n r 



+o 

where r c is the zeroth-order distance of closest approach defined as 
r r = r B |iV R x nJ. 



(49) 



(50) 



A similar reasoning applied to (|38fe[) yields the expression of the triple l r (N e , x B ) 
characterizing the direction of r s (iV e , x B ) at x B as a function of N e and x B . 

Proposition 6 Under the assumptions of proposition the triple yielding the 
direction of propagation of a quasi- Minkowskian light ray r s (7V e , x B ) at point x B is 
given by 

T r (N e ,x B ) = -N e 



m 
r r 



(l+7)|JV e x n B \ + 



n\N e xn B \ 2 - (1+^(1 + N e .n B ) 



>N P 



lib I . , . 1 1 I \ 

+ — { (1 + + N e .n B ) + —{k 



7t — arccos(AT e .n B ) 



IN- x n B \N e .n f 



(l + N e .n E 
{ +7j N P xn B 



(51) 



where P e is the unit 3-vector orthogonal to N e defined as 
„ N e xn B 

x N P . 



(52) 



\N e x n E : 

Of course, the unit vector N e and the impact parameter b may be taken as the 
boundary values which determine a quasi-Minkowskian light ray emitted at infinity. 
Then the notation T s (N e ,b) instead of T s (N e ,x B ) will be used. The following 
proposition is straightforwardly deduced from (141 

Proposition 7 Let T s (N e ,b) be a quasi-Minkowskian light ray emitted at infinity in 
the direction N e with an impact parameter b. At any point x B of T s (N e ,b) fulfilling 
condition Cjfl), the direction triple is given by 



l r (N e ,x B ) = -N e - 



m\N„ x n, 



1+7 + 



kto|7V„ x n F 



ill i , . , _ . ruiiL 

+-rl(l+l)(l + N e .n B ) + — 



7r — arccos(AT e .n B ) 



IJVfi x n B \N„ .Ti f 



O 



b 3 



(53) 



We note that using an expansion in powers of m/b instead of m/r c cancels the 
'enhanced' term — (1 +j) 2 (m/r c ) 2 (l + N e .n B ) 2 /\N e x n B \ appearing in the right-hand 
side of (fBTj) . We shall discuss this question in the next section. 
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7. Bending of light rays coming from infinity 

In order to define the gravitational deflection of a light ray relative to a static observer 
S B staying at a given point x B , let us consider first two quasi-Minkowskian light 
rays T s and r' s arriving at the same instant at x B . We suppose that T s and T' s 
are emitted at infinity in initial directions N e and N' e , respectively. Furthermore, 
we assume that a condition like (|48|) is met by each of these rays. If gravity was 
lacking, the angle between T s and 1^ as measured by S B would coincide with the 
Euclidean angle 8 £ (N e ,N' e ) formed by N e and N^. The value of this angle may be 
considered as an intrinsic quantity since N e and could be determined (at least in 
principle) by measurements performed by static observers located at spatial infinity. 
As a consequence, the gravitational contribution to the angular separation between 
T s and T' s as measured by S B may be characterized in an intrinsic way by a definition 
as follows. 

Definition 3 Suppose that a static observer at point x B is receiving two quasi- 
Minkowskian light rays T s and T' s emitted at infinity in the directions defined by the 
unit vectors N e and N^,, respectively. The gravitational contribution to the angular 
separation between T s and T' s as measured by this observer is by definition the angular 
quantity 

6i gr) (T s ,r' s ) =S B (T s ,T' s )-S E (N e ,Ni). (54) 

Taking (|22p into account, it is immediately seen that 6 B 9r \r s , T' s ) may be 
determined by the relation 

6 { B 9r) (T S: T' s ) = arccos [A(r B )B(r B ) T r {N e , x B )d r (N' e , x B )] - arccos(N e .N' e ). (55) 

Definition [3] and equation (1541) enable to obtain a coordinate-independent 
characterization of the gravitational deflection of a quasi-Minkowskian light ray coming 
from infinity. There exists one and only one null radial geodesic emitted at infinity and 
arriving at x B at the instant of reception of T s . Pointing out the quasi-Minkowskian 
nature of this ingoing null radial geodesic, we may use definition [3] in order to propose 
the new definition below. 

Definition 4 Let T s be a quasi-Minkowskian light ray emitted at infinity in a given 
direction N e and observed at a point x B fulfilling condition [48\j - The gravitational 
deflection of this ray relative to a static observer at x B is by definition the angular 
quantity 

A XB (r s ) = s { r } (T s ,T; d ), (56) 

where T~ d is the ingoing radial null geodesic path arriving at x B at the instant of 
reception of T s . 

The sign occurring in the right-hand-side of (|56[) is chosen in order to have a 
positive quantity for the light deflection. The ingoing null radial geodesic T~ d is 
emitted at infinity in the direction given by the unit vector — n B . So the direction 
triple of this geodesic at x B may be denoted by l r (—n B ,x B ) according to the notation 
already used in section [5] A closed-form, exact expression of l r (— n B ,x B ) is easily 
derived from the geodesic equations. One gets 
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Taking ([571) into account, it is easily seen that ([2~Tj) implies 5 s (r s ,r~ d ) = 
5 £ {l r {N e , x B ), n B ). Hence an equation as follows, directly inferred from equations 
HMH and (JSU): 

A Xfl (r s ) = S s (-N e ,n B ) - 5 E (l r (N e ,x B ),n B ). (58) 
Since l r (N e ,x B ), n B and N e are coplanar triples, one has 

8 £ (-N e ,n B ) = 5 e (-N e ,L r (Ne,x B )) + 5 £ (l r (N e ,x B ),n B ). (59) 
Substituting for 5 £ (— N e , n B ) from (I5T)1) into ([55| . we obtain the following proposition. 

Proposition 8 The gravitational deflection of a quasi- Minkowskian light ray 
T s (N e ,x B ) relative to a static observer at x B is equal to the Euclidean angle between 
—N e and the direction triple of the ray at x B : 

A Xb (T s ) = S e (-N e ,l r (N e ,x B )). (60) 

In any practical case encountered in high-accuracy astrometric projects, the 
gravitational deflection is a small positive angle. As a consequence, (|6T)]) is equivalent 
to 

A XB (r s ) = arcsin A(r B )B{r B )\T_ r {N e ,x B ) x 2V e |) • (61) 

Equations (f54"|) - (f6"Tj) are exact relations. Within the post-post-Minkowskian 
approximation, however, the terms of order G 3 may be neglected. So the expression 
of the gravitational deflection reduces to 

A XB (r s ) = y/A(r B )B(r B ) \l r (N e ,x B ) xN e \ + . (62) 

Substituting for l r (N e ,x B ) from (|55|) into (|6"2"j) . and then replacing m by GM/c 2 
in order to facilitate precise numerical estimates, we get a proposition as follows. 

Proposition 9 Let T s be a quasi- Minkowskian light ray emitted at infinity in a given 
direction N e and observed at a point x B fulfilling condition \4<ty - The gravitational 
deflection of T s relative to a static observer at x B is given by 

(1 +7)GM, 

~c 2 h~ 



A Xb (T s ) = V (1 + cos</> B ) 



G 2 M< 



c 



^ sin20 B I — (1 + 7)^(1 + cos 4> B ) sin ( 



2 



G 3 Af 3 



(63) 

where b is the impact parameter of the ray determined by equation J^ff| ) and 4> B is the 
Euclidean angle formed by N e and n B , namely the angle defined by 

N e .n B = cos cf) B , < 4> B < it. (64) 



The term of order G in (1631) is currently used in VLBI astrometry. The 
contribution of order G 2 may amount to about 5.5 microarcseconds (/ias) for a ray 
coming from infinity and arriving tangentially at the surface of the Sun. So the terms 
in G 2 will be indispensable in a foreseeable future in order to perform highly precise 
tests of general relativity with laser rays passing near the Sun (LATOR or ASTROD 
projects, e.g.). 
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The right-hand side of equation (|63[) yields a coordinate-independent expression 
of the deflection of light which involves only two quantities, each of them having a 
perfectly clear geometrical meaning. Indeed, 6 is a length which could be measured 
by a static observer at infinity, and 4> B is the angle formed by two directions defined 
in the rest space at infinity. Nevertheless, the intrinsic nature of Ax B (r s ) does not 
automatically prevent some problems from raising in the discussion of the different 
contributions. To see that, let us replace b by its expression yielded by equation (|49l) . 
Then Ax B (T s ) reads 

(1 + 7 )(7M 



Ax B (r s 



-(1 



G 2 M 2 



-(1 +7) 2 (1 + cos0 B ) 



— sin2<?!> f 
2 



COS ( 



sm < 



O 



( G 3 M 3 
V c 6 r 3 



(65) 



where r c — r B sin 4> B . Let us suppose that we apply (|65[) to a ray deflected by a 
massive body of the Solar System and observed by Gaia, which occupies an orbit 
about 1.5 million kilometres from the Earth around the Lagrange point L2. Denoting 
by r P the value of r at the periapsis of the ray, we have r c « r P and sin0 B rs r P /r B . 
Then the term in 1/ sin <p B appearing on the right-hand side of equation (|65|l seemingly 
generates a post-post-Minkowskian contribution 

Ax ^n,) RJ _ (1+7) 2 (1 + cos ^ ) 2" Tb (66) 

crrp r P 

which is equivalent to the so-called 'enhanced' post-post-Newtonian term given by 
equation (94) of [TO]. The magnitude of this contribution is totally negligible when 
the deflecting body is the Sun. Indeed, the fact that in this case r B « 1 AU and 
r P > 0.7 AU for any light ray observed by Gaia implies | Ax B '°" h) \sun < 6.7 x 10 _4 /Ltas. 
Nevertheless, the right-hand side of (|6"6"|) may be comparatively large, e.g., for a 
light ray grazing Jupiter; one has indeed in this case 4AU < r B < 6AU whereas 
r P w 7.14 x 10 4 km. Then 10.7^as < | Ax B '" nh> [jupiter < 16/ias, which means 
that the value of the effect may be slightly greater than the standard error about 
10 /zas expected for the brightest stars during the Gaia mission [21] . This effect is 
illusory, however, since r c differs from r P by a term of order G which compensates the 
contribution (|66[) , as it is shown by the absence of any term in 1/ sin0 B in the right- 
hand side of equation (|6"5)) . So the 'enhanced' term Ax B ' c " h) must be considered as a 
fictitious contribution due to the fact that neither r c nor r P are intrinsically defined 
parameters. This analysis confirms the conclusions drawn in [10] and [22] , 



8. Total deflection of light 

It is now easy to carry out the calculation of the total deflection Ax of a quasi- 
Minkowskian light ray T s emitted at infinity with an impact parameter b 3> r/j and 
returning to infinity after passing near the central body. Since condition (|48p cannot 
be fulfilled for any point x B on T s , we shall restrict our attention on the portion of 
T s between the emission point at infinity and the periapsis of the ray, say x P . Owing 
to the symmetry, we can put 

A X = 2A Xp (T s ), (67) 
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where Ax P (r s ) denotes the gravitational deflection of the light ray T s (N e . x P ) relative 
to a static observer at x P , N e being the emission direction of T s at infinity. Since 
l r (N e ,x P ).n P — 0, we have 

5 £ {Lr(Ne,x P ),n P )=n/2. (68) 

Substituting for 5 e (l r (N e ,x P ),n P ) from (|68|) into ([58)) written at x P , noting that 
5 £ {— iV e , n P ) — 7r — P , where <fi P is the Euclidean angle between N e and n P , and 
then taking (|67j) into account, we find 

m 

Replacing now A% B (r s ) by Axp(r s ) and <ft B by 4> P in ([63| . and then taking into 
account ([67]) and (|69|) . we get 

G 2 M 2 



Ax 



2(1 + 7 )GM 



c 2 6 



1 + sin 



-2(1 + 7) 2 



Ax 
2 

1 + sin 



c 4 6 2 



k (tt + Ax + sin Ax) 



A X 



/ G 3 M 3 
\ c 6 6 3 



(70) 



Finally, solving (|T0[) for Ax, and then neglecting all terms of order higher than m 2 /b 2 , 
we are led to the proposition which follows. 

Proposition 10 The total gravitational deflection Ax of a quasi- Minkowskian light 
ray starting from infinity with an impact parameter 6 3> r^ is given by 

2(1 + 7 )GM kttG 2 A/ 2 „ / G 3 Af 3 " 



A X = 



~ 2 b 



» 4 b 2 



O 



c% 3 



(71) 



Equation (fTTj) agrees with the results classically yielded by integrating the null 
geodesic equations (see, e.g., [TU] and references therein). 



9. Conclusion 



In this paper, the direction of light propagation is obtained at the post-post- 
Minkowskian level for a three-parameter family of static, spherically symmetric 
spacetimes without integrating the geodesic equations. It must be emphasized, 
however, that the explicit expansions obtained in sections [S][5] only hold for the null 
geodesic paths that we have called the quasi-Minkowskian light rays. It may be noted 
that this limitation is also encountered in the previous works treating the problem by 
integrating the differential equations of null geodesies. 

Our method based on the time transfer function is perfectly adapted to the generic 
case where the emitter and the receiver are both located at a finite distance from the 
source. The central result of the paper is constituted by equations (|38a|) and (|38b|) 
yielding explicit expressions for the direction triples of a quasi-Minkowskian light ray 
at its emission and reception points. Supplemented by the formula (|36[) yielding the 
impact parameter b of the ray, these equations straightforwardly lead to equations 
(|41al) and (141 b[) giving the expansions of the direction triples in powers of m/b. 

The expansions obtained in the generic case are easily extended to a light ray 
emitted at infinity, the receiver being an observer staying at a finite distance from 
the centre. An intrinsic characterization of the gravitational deflection angle of such 
a ray observed by a static observer is formulated in definition |3l The well-known 
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expression of the deflection angle currently used in VLBI astrometry is extended to 
order G 2 by equation (|(53")) . A recent discussion relative to an apparent 'enhanced' 
post-post-Minkowskian term is confirmed. Finally, the well-known formula giving the 
total light deflection angle up to order G 2 is recovered. 
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Appendix A. Conservation of the angular momentum 



In order to prove that the angular momentum L defined by equation (lll[) is conserved 
along any geodesic path of the metric ([1} , it may be assumed without lack of generality 
that £ is an affine parameter. The geodesic paths is then governed by the system of 
Euler-Lagrange equations 



dc, 



dlj 



A'(r)(l°f 



where 



B 2 {r) 

dx^/dC. Equation (|A.1|) implies 

dx° 



B'(r) kl 

Okll I 



X 

r 



l = A(r)- 



d( 



E, 



(A.l) 
(A.2) 

(A.3) 



where E = const. Excluding now the spacelike geodesies located in an hypersurface 
x° = const, we can suppose that E ^ 0. So dividing (|A.2[) by Iq an d taking (IA.3|) into 
account, it is easily seen that the geodesic path satisfies an equation as follows 



dl 



E 
~2 



A\r) 
A 2 (r) 



B'(r) (l 



(A.4) 



where I is defined by equation ©. Equations dl]) and (|A.4|) imply the conservation 
law d(—x x l)/dQ = 0. 
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